Extended Einstein relations with a complex effective temperature in a one 

dimensional driven lattice gas 
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We carried out numerical experiments on a one-dimensional driven lattice gas to elucidate the 
statistical properties of steady states far from equilibrium. By measuring the bulk density diffusion 
constant D, the conductivity a, the intensity of density fluctuations \y we confirm that the Einstein 
relation D\ — aT, which is valid in the linear response regime about equilibrium, does not hold 
in such steady states. Here, T is the environment temperature, and the Boltzmann constant is 
set to unity. Recalling that the Einstein relation provided the first step in the construction of 
linear response theory, we attempt to extend it to a generalized form valid in steady states far from 
equilibrium. In order to obtain new relations among measurable quantities, we define a complex 
effective temperature B — i$ from studying the static response of the system to a slowly varying 
potential in space. Replacing T in the Einstein relation by the real part of the effective temperature 
0, we numerically confirm that the relation D\ = o~Q holds in the nonequilibrium steady states 
far from equilibrium that we study. In addition to this extended form, we find the new relation 
(L/2n)cx = cr$, where c represents the propagation velocity of density fluctuations. 

PACS numbers: 05.40.-a, 66.30.Dn, 02.50.Ey 



I. INTRODUCTION 

In 1905, Einstein derived the relation D = \±T between 
the diffusion constant D for non-interacting Brownian 
particles suspended in an equilibrium fluid of tempera- 
ture T and mobility /j with respect to an external force 
0. (Here, the Boltzmann constant is set to unity.) This 
relation is now called the 'Einstein relation'. Through 
the development of linear response theory, the original 
Einstein relation was extended to the form 



(1) 



which is valid for interacting many-body systems Q. 
Here, D should be considered as the bulk density dif- 
fusion constant, x is the intensity of density fluctuations, 
and a is the conductivity. 

The Einstein relation is a universal relation in the lin- 
ear response regime about equilibrium. It represents the 
first step in the construction of nonequilibrium statisti- 
cal mechanics in the linear response regime. With this 
in mind, if the Einstein relation can be extended to some 
form valid for nonequilibrium steady states (NESS) far 
from equilibrium, realizing an understanding of this form 
will be an important step in the construction of statistical 
mechanical theory of NESS. 

With this motivation, because it is simpler to study 
one-particle NESS systems than many-body systems, we 
previously analyzed a one-particle Langevin equation 
with a tilted periodic potential Q- In that study, it 
was found that there exists an extended Einstein rela- 
tion with an effective temperature that can be measured 
independently of the diffusion constant and the differen- 
tial mobility. Then, it is natural to ask whether that 
result for a Brownian particle in NESS can be extended 
to many-body interacting systems. In the present paper, 
we address this question by studying a one-dimensional 



driven lattice gas (DLG) 0, |H El as an example of a 
many-body interacting system. 



II. MODEL 

Let r\ x be an occupation variable defined on each site 
a; of a one-dimensional periodic lattice x = 1, 2, • • • , L. 
The variable r\ x takes the value 1 if the x-th site is oc- 
cupied by a particle and if it is empty. In order to de- 
scribe the model simply, we introduce the variable j/l+i, 
which is assumed to be equal to r\\. The array of occu- 
pation variables (771, • • • ,rjr,) is denoted by rj and called 
the "configuration" . 

The time evolution of 77 is described by the follow- 
ing rule: At each time step, randomly choose a near- 
est neighbor pair x and x + 1, where 1 < x < L, 
and exchange the values of rj x and ffc+i with probabil- 



ity c(x; rj) given by c(x; rf) 



ae 



where rf is 



the configuration obtained from r\ through this exchange, 
P = 1/T is the inverse temperature of the environment, 
and Q(r) — > rf) represents the heat absorbed from the 
heat bath as a result of change of configuration rj —> i] x 
in one Monte Carlo step. The constant a is determined 
so that c(x;r]) < 1. The particle number N = Y] r) x 
is conserved throughout the time evolution, and the av- 
eraged density p — N/L is a parameter of the model. 
In this paper, we study the DLG for which 



Q( V - rf) = H (rj') - H ( v ) - E^Uv 



Tj 



(2) 



x = l 



where E is a constant driving force, Ho(r)) is the equilib- 
rium Hamiltonian H (rj) = - J2 x =i VxVx+i, and jx(v ~ * 
rf) is the net number of particles hopping from x to x + 1 
for the configuration change rj — > rj' in one Monte Carlo 
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FIG. 1: j(t) is plotted with error bars in the case p = 0.5, 
L = 64 and \Vi\/VL = 0.1 (See JTJ3J). The slope of the 
fitted line (the dotted line) in the early stage (t < 150 mcs) 
corresponds to D(2ty/ L) 2 . The slope of j(t) in the late stage 
(t > 150 mcs) differs from D(2iv/L) 2 . 



FIG. 2: D as a function of L for the case p = 19/64 (solid 
circles) and p = 0.5 (open circles). The dotted line and the 
solid line represent 0.042\/L + 0.20 and 0.068^ + 0.0413, 
respectively. 



step. In the argument below, ( ) E represents the statis- 
tical average in a NESS with driving force E. 

Throughout this paper, we fix /3 = 0.5, and the numer- 
ical parameter a introduced in the exchange probability 
c(x;r]) is set to exp(— 3/3). This value is chosen so that 
c(x,v) < 1 in the case E < 5. As a typical example of 
a NESS far from equilibrium, we study this model with 
E = 4, and we consider the cases p = 19/64 and 0.5, 
with L = 64, 128 and 256. We remark that because we 
attempted to evaluate precisely measurable quantities in- 
troduced in the next section, it was hard to obtain data 
in the case L > 256. 



III. EINSTEIN RELATION 

Because we have little knowledge regarding relations 
characterizing NESS of DLG, let us start by investigating 
the validity of the Einstein relation Q by defining the 
quantities a, x an d D for this model. First, using the 
averaged current Je, which is defined by 



J e = j (J^jxiv -> V') 



the conductivity a is expressed as 

dJ E 



dE 



(3) 



(4) 



Second, the intensity of large scale density fluctuations 
is defined by 



where 



x(L) = (\f ll \ 2 ) s E , 



fj n (t) = ^=X>(£f)e <2pa! - 

V X — l 



(5) 



(6) 



Third, the density diffusion constant D is defined as the 
coefficient of the diffusion term of the evolution equation 
describing the averaged behaviors of the density field. 
As an example of such behaviors, we consider the time 
evolution of the density field after a slowly varying weak 
potential V x is turned on at t = in the case that the 
system is in a steady state. Here V x is given by 



V x = -^=Vie~ 



c.c, 



(7) 



where c.c denotes the complex conjugate, and V\ is a 
complex number satisfying \V\\/\/L <C 1. 

The averaged behavior in this case is represented ap- 
proximately by the time evolution of a quantity (^i(i))^, 

where ( ) E denotes the statistical average under the 
condition that the equilibrium Hamiltonian becomes 

Hv(v) — Hoiv) + Ylx=i^^> when t > 0. We define 
the response function against V x as 



|Vi|—o Vi 



(8) 



In order to see the behavior of R(t), in Fig. ^ we plot 



7(f) = log 



V 



(9) 



Note that we have checked that 7(i) takes the same value 
for the cases |Vi|/VX = 0.05, 0.1, 0.15 and 0.2. It implies 
that we can estimate the value of the right hand of © 
within a numerical accuracy. It is seen in Fig. ^that j(t) 
has the linear slope in the early time regime. Therefore, 
using this linear slope, D is estimated from the form 



|.R(oo) - R(t)\ = const.e- 15 ^) 2 *, 



(10) 
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FIG. 3: The profiles of {rix) s E v under the influence of V x — 
Asm(2Tvx/L) in the case p = 19/64 and L — 64. A = 0.1, 
0.2, 0.3 and 0.4 correspond to the circle, square, triangle and 
plus, respectively. 



in the early time regime. (We remark that we cannot 
deny a possibility that the linear slope shown in the early 
time regime might be caused by finite size effects. This 
is an unresolved problem.) 

Using the above definitions of er, x an d D, we con- 
firm the validity of numerically for the equilibrium 
case, E = 0. With regard to other relations in the linear 
response theory, we can derive 0) from the Green-Kubo 
relation and the fluctuation response relation, which have 
been proved for DLG models y| ■ Contrastingly, in NESS 
for our DLG, we find that the two sides of differs 
greatly. Specifically, in the case p — 19/64, we have con- 
firmed that tr(L) and x{L) converge to er* = 0.032 and 
X* = 0.26 in the limit L — oo, while D has the anoma- 
lous L dependence D(L) cx yZ, as shown in Fig. |2] 
This anomalous behavior is essentially the same as that 
reported in Ref. Also, in the case p — 0.5, is 

not valid because we have confirmed that cr* = 0.037 and 
X* = 0.32 in the limit L — > oo, while D(L) cx \L (See 
Fig. EJ). 

Now we ask whether it is possible to obtain an exten- 
sion of the Einstein relation valid for NESS of interacting 
many-body systems by defining an appropriate measur- 
able effective temperature. This question arises natu- 
rally, because we have found that such an extended Ein- 
stein relation does hold for NESS of one-particle systems. 
In this relation, the temperature of the environment, T, 
is replaced by an effective temperature, which can be 
measured independently of the other quantities appear- 
ing in the Einstein relation Q . By analogy to that case, 
we seek to define an effective temperature which when 
used in place of T in Q provides an extended Einstein 
relation that is valid for NESS in interacting many-body 
systems. 
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FIG. 4: The profiles of {r)x) s E v under the influence of V x = 
Asin(27ra;/L) in the case p = 0.5 and L = 64. A = 0.1, 0.2, 
0.3 and 0.4 correspond to the circle, square, triangle and plus, 
respectively. 

IV. EXTENDED EINSTEIN RELATIONS 

In a previous study, we showed that the effective tem- 
perature appearing in the Einstein relation valid for 
NESS of a Brownian particle can be measured by adding 
a slowly varying potential to the system In a simi- 
lar way, we wish to determine an effective temperature 
by adding to the system a slowly varying potential given 
by 0. Due to the influence of this potential, we have 
(fii) s E v 7^ 0, where {}% v denotes the steady state aver- 
age in the system with driving force E and potential V x . 
Recalling that the relation 

{flx) s E = ,v = -Vi^ + o{\V l \), (11) 

holds in the equilibrium case, we define the effective tem- 
perature B\ in NESS as 

- lim ^ = m. (12) 
| ft |—o Vi 0i 

Note that in NESS, the left-hand side of |JT2J| is a com- 
plex number in general. To see this concretely, we 
show that (r) x ) E v under the influence of the potential 
V x = A sin(27rx/L) has a component of cos(2ttx/L) (see 
Fig. EJ) • Expressing explicitly that 9\ is a complex num- 
ber, we write 

§x = e-i<f>, (13) 

where O and $ are real numbers. In Fig. [SJ and <& 
are displayed as functions of L in the case p — 19/64 
and p = 0.5. It is seen that we have the anomalous L 
dependences @(L) cx s/L and <&(£) cx L. We remark that 
$ = in the case p = 0.5 (See Fig. EJ. 

Noting that from the graphs in Fig. Inland Fig. |SJ D(L) 
and Q(L) are both approximately proportional to y/~L, we 




FIG. 5: O and $ as functions of L in the cases p — 19/64 (solid 
symbols) and p = 0.5 (open symbols). Q(L) and corre- 
spond to the circle and square, respectively. The solid line, the 
dotted line and the dash-dotted line represents 0.36\/L + 1.4, 
0.21L and 0.58^1 + 0.26. Note that T = 2. 



FIG. 6: D% as a function of <j6. The cases p — 19/64 and 
p = 0.5 correspond to the solid symbols and the open symbols, 
respectively. The cases L = 64, 128 and 256 correspond to the 
square, triangle, and circle, respectively. The straight line 
represents D\ = crQ. 



investigate whether the Einstein relation can be extended 
to the NESS by simply replacing T in Q by 8(L). In 
Fig. El we plot D(L)x(L) as a function of a(L)Q(L) 
for the cases p — 19/64 and 0.5, with L = 64,128 and 
256. The graph clearly shows that the extended Einstein 
relation 



D(L) X (L) = a(L)e(L) 



(14) 



holds within the precision of the numerical computations 
0. 



V. FAILURE OF EXTENDED FLUCTUATION 
RESPONSE RELATION 

There are universal relations in the linear response 
regime about equilibrium, such as the Einstein relation, 
the fluctuation-response relation (FRR) and the Green- 
Kubo relation. One may expect that other relations can 
be extended in the similar way as did for the Einstein re- 
lation l|14|) in NESS far from equilibrium. In this section, 
we investigate an extension of the FRR. 

The FRR connects the time dependent response func- 
tion and the time correlation function. In order to define 
the response function, we consider the situation that the 
perturbation Q is applied to the system when t = 0. 
The response function R{t) in this situation is defined by 
In the equilibrium case, we can prove the FRR 



1 



R(t) = -C(t), 



with the time correlation function 



C(t) = ^(\m(t)-fj 1 (Q)\ 2 Y 1 



(15) 



(16) 



Because the Einstein relation was extended with the 
replacement of T by 0\ , we expect that Ijl5|l can be ex- 
tended in the form 



R(t) 



1 



C{t). 



(17) 



Specifically, we numerically investigate the validity of 
(I17|l in the case p = 0.5 and L = 64. Because $ = 
in the case p = 0.5, 1)17(1 is rewritten as R(t) = Q~ 1 C(t). 
In Fig. [7| we plot R(t) as a function of C(t) during an 
interval < t < 800 mcs. It is seen that R(t) devi- 
ates from &~ 1 C(t) substantially in an early time regime 
(0 < t < 150 mcs) where the diffusion constant D was 
defined (See ltTU|) ). This implies that ((T7jl is not valid. 
Note that R(oo) = 8~ 1 C(oo) holds by the very definition 
of 6>i given by {Bl) . 

We remark that a generalized version of FRR of the 
nonlinear Langevin dynamics without the detailed bal- 
ance condition was presented by several authors |17| . 
However, in their formalism, a generalized force which 
defines a response function is chosen so that the sim- 
ple form of FRR holds, and it seems that their response 
function cannot be measured directly in experiments. In 
contrast to the previous studies, we consider a relation 
between C(t) and R(t) which are measured directly. 



VI. INTERPRETATION 

In sections Hvl andlVl we find that the Einstein relation 
(|14|l is extended with the diffusion constant D which is 
measured from R(t) in the early time regime and the 
effective temperature O, while the extended FRR <|17|) is 
not valid unfortunately. In paragraphs below, we shall 
provide the phenomenological understanding of validity 
of {Tl} and invalidity of ((T7jl . 
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FIG. 7: R(t) as a function of C(t) in the case p = 0.5 and 
L = 64 during an interval < t < 800 mcs. The slope of the 
thin line represents O -1 . 



The extended Einstein relation 1)14(1 can be explained 
from the following argument. Let J(p, E) be the av- 
eraged particle current in a system with density p and 
driving force E. When the slowly varying potential Q 
is applied to this system, first the current should become 
J((Vx)e v ' E ~ dV x /dx), and as a result a diffusive cur- 
rent D ((jfa+i — Vx))e v should appear so as to make the 
total current homogeneous in x; that is, it is reasonable 
to assume the quantity 



J (Vx) E y,E 



Vx)) 



E,V 



(18) 



is independent of x in a steady state. This can be rewrit- 
ten as 

/as b ,. 27T- 2tt B (ivil (m) E ,v \ 
c K~ni) E ,v+ l(J -^ 1 ~L~ E - y = 1 T - ' — T~ I ' 

(19) 

where c = dJ{p, E)/dp represents the propagation veloc- 
ity of density fluctuations. From 1|12[) and l(19|l . we obtain 



a0 x = Dx-i—xc- 

Z7T 



(20) 



This leads to the extended Einstein relation (|14fl . 

Here, we note that a microscopic derivation of l|18|l is 
quite difficult. In particular, we have not been able to 
show that D defined by (|1L)|> and D in (|18fl are identical. 
In relation to this point, we note that a different value 
of D is obtained if we consider the time regime t > 150 
mcs in Fig. ^ Using this value of D, the relation 114|) 
cannot be confirmed numerically. This result suggests 
that memory effects in a single mode description cannot 
be neglected in the late stage of the relaxation process 
after the potential is turned on. Because the memory 
effects arise from a nonlinear coupling between different 



FIG. 8: B(t)t as a function of t in the case p — 0.5 and L = 
128. The guided line (the dashed line) represents B(t)t oc t A ^ . 

modes, the simple assumption l(18JI in NESS is far from 
the trivial. 

The non-trivial nature of the extended Einstein re- 
lation l|14f) can be seen through the failure of the ex- 
tended FRR 1)17(1 . For example, because we cannot ob- 
serve the exponential decaying regime of C(oo) — C(t) 
even in the early time regime unlike the behavior of 
R(oo) — R(t) (See Fig. HJ, we cannot define the dif- 
fusion constant D properly from the assumption that 

log \C(oo) - C(t)\ = -D (^) 2 t + const., though in the 
equilibrium case, the same value of D is obtained whether 
it is defined from log |C(oo) - C(t)\ or log \R(oo) - R(t)\ 
because of the FRR lfT5]l. 

The non-exponential decaying behavior of C (t) may be 
related to a long time tail of the time correlation func- 
tion of current fluctuations. Indeed, using the continuity 
equation, C(t) is expressed by 



\ i, - 1. 1 - 1 

(21) 

It is expected that the behavior of C(t) in a shorter time 
regime than the diffusion time D^ 1 (2tt/L)^ 2 has com- 
mon features with that of 



B(t) = \ 



Li 



EEkw*- 1 )-^*))-- 7 *] 



k=l x=l 



E 

(22) 

which is displayed in Fig. [HI in the case p = 0.5 and L = 
128. The graph is consistent with the relation B{t)t oc 
t 4 / 3 which was derived from the theoretical analysis of 
fluctuating hydrodynamics 

We remark that these observations on R(t) and C(t) 
cannot be explained from the assumption that the time 
evolution of density field is described by a linear Langevin 
equation of rji (t) , because under this assumption both the 
relations l|14|) and i|17|) are derived though l|17|) was found 
not to valid. This implies that the dynamical aspect of 
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FIG. 9: k~ 1 cx as a function of cr$, where k — 2-k/L, for the 
case p = 19/64. Note that for the case p — 0.5, c = $ = 0. 
The cases L = 64, 128, 192 and 256 correspond to the square, 
triangle, diamond and circle, respectively. The straight line 
represents (L/2n)cx = cr$. The error bars represents the 
standard deviations of the statistical samples. In addition to 
the statistical error bars, there exists a fitting error which is 
less than 5 percent of the quantities. 



the system is too complicated. It may be surprising that 
the extended Einstein relation holds with the appropriate 
choice of D. 

Finally, from (|20[1 . we find a new relation 



^c(L) X (L) = *(L)$(L). 

Z7T 



(23) 



We investigated the validity of (|23[1 in numerical exper- 
iments. We numerically measured c from the oscillatory 
behavior of the time correlation function i|16|) . In one 
example of this analysis, we considered all ti satisfying 
C{ti) — C(oo) = 0, where < ti. Then, c can be 
estimated from the relation 



2tt , . 
c— [ti+i - ti) 



(24) 



We have confirmed that c converges to c* = 0.17 in the 
case p — 19/64 in the limit L — > oo. In Fig. |§1 we 
plot (L/2tt)c(L)x(L) as a function of a(L)$>(L) for the 
cases p = 19/64, with L = 64, 128 and 256. (Note that 
c = $ = for the case p — 0.5.) This graph clearly shows 
that the extended Einstein relation H23|) holds within the 
precision of the numerical computations. 



VII. THERMODYNAMIC CONSIDERATION 

Here, we devote space to thermodynamic consideration 
related to our results. We first study the equilibrium case 
E = by considering a situation that a perturbation 
potential is applied to the equilibrium system. When we 
choose a slowly varying potential, we can regard V x and 



(Vx) S e =0 v as smooth quantities V{x) and p(x). By use of 
the chemical potential /x(p, T) of the equilibrium system 
(without the potential V(x)), a density profile under the 
influence of V(x) is determined by a relation that 



/i(p(x), T) + V(x) = const, in x, 



which is equivalent to 

dp\ dp dV g 
dp J T dx dx 

Note that the fluctuation relation 

-l 



dp y T 



X 
T 



(25) 



(26) 



(27) 



holds in the equilibrium state. As another relation in this 
situation, we have the total zero current condition 



dp -( dV 
ax \ ax 



0. 



(28) 



By expanding the second term of (|28|) in dV/dx, we ob- 
tain 



dx dx 



(29) 



A combination of (HSJ , JSJJl and leads to the Einstein 
relation In the equilibrium case, we can show just 
the opposite that a combination of the Einstein relation 
and the relation (|29|l derives the fluctuation relation 
(|2Tjl through the relation 1(1251) . 

In connection with the above mentioned thermody- 
namic consideration, we introduce our previous study 
on a two-dimensional DLG ^4|- First, we remark that 
an appropriate definition of the chemical potential for a 
NESS is not known, unless we use the local equilibrium 
assumption. Because in NESS far from equilibrium, the 
local equilibrium assumption is not valid, we attempted 
to define a chemical potential without this assumption. 
Focusing on a density response in a direction perpendic- 
ular to the driving force, we defined a nonequilibrium 
chemical potential by a similar relation with (|25|l and 
constructed an extended thermodynamic function from 
measurable quantities. As a non-trivial prediction, we 
showed numerically that the rate function of density fluc- 
tuation was determined by this extended thermodynamic 
function. A theoretical framework behind these observa- 
tions has been uncovered We remark that the Ein- 
stein relation Q was also observed in a direction perpen- 
dicular to the driving force |l6j| and that the argument 
in the previous paragraph can be applied to this case. 

In the present paper, we have reported that the ex- 
tended Einstein relation (|14|> and the new relation 123|) 
hold in NESS, and a combination of both the relations 
and (|19|) that corresponds to the Fourier transform of 
(|29|l leads to the definition of the effective temperature 
l|12p. Because in the case E = 0, the relation (|11|) is 
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FIG. 10: In the case H (rf) = (ASEP), j(t) is plotted with 
error bars in the case p = 0.5, L — 64 and |Vi|/\/L = 0.1 (see 
llOl l'l. The slope ol the fitted line (thin line) in the early stage 
(t < 150 mcs) corresponds to D(2ir/ L) 2 . 

derived from (|26[l and l|27l) . one might expect that there 
exists a nonequilibrium chemical potential even for the 
present NESS. For example, as a trial, one might propose 
a nonequilibrium chemical potential jli(pi,9i) as 

£iG3i,0i)+Vl =0. (30) 

However, |5DJ and (JT^Jl lead to 

{-Y-f- < 3i > 

Because the relation (|31|) is hard to be interpreted as an 
extended form of the fluctuation relation Q27p. the defini- 
tion (|30l) seems not to be useful. Taking into considera- 
tion that the relation (|30l) is just one trial, there might be 
a plausible definition of a chemical potential for NESS, 
but we have not found it yet. The difficulty that does 
not arise in the equilibrium case lies in the wavenumber 
dependence of D, as is expected from the system size 
dependence of D shown in Fig. [5J This wavenumber 
dependence is not observed in the equilibrium case and 
prevents us from pushing ahead further thermodynamic 
consideration. 



VIII. DISCUSSION 

Equations l|14|) and (|23|l constitute the main results of 
this paper. In the following, we give five important re- 
marks that are peripherally related to these main results. 

First, we give an argument from which to understand a 
physical meaning of 9\ as the effective temperature. Ef- 
fective temperatures have been studied mainly in glassy 
systems jy, El Ell [K| . The simplest characterization 
of temperature may be provided by the equi-partition. 
The law of equi-partition is regarded as one example 



of an energy distribution law that can be expressed by 
the canonical distribution. Out of equilibrium, the equi- 
partition law and the canonical distribution do not hold 
in general. Indeed, the steady distribution p s E (rj) dif- 
fers greatly from the canonical one. However, as demon- 
strated in Refs. 0, IT?! ], by constructing measurement 
methods to extract behavior at some large time and/or 
length scales, the equi-partition law and the canonical 
distribution are recovered at these scales. 

Hypothesizing a similar situation for the NESS we 
study, we assume a simple form for the modification of 
the steady state measure under the condition that the 
potential varying slowly in space is added to the system. 
The form we consider is 

P s e.v(v) = ^(^e-W-^^+c.c) (32) 

ZjV 

for | VI | <C 1, where Zy is a normalization constant. From 
(|32|l . we can derive (|12l> . which is the definition of 0\. 
We believe that the simple form Ij32(l . which describes 
the response of the steady state measure to this slowly 
varying potential, is a good starting point to investigate 
NESS. 

Second, we remark that our main result also holds for 
the case Hq(t]) — 0, which is an example of asymmetric 
simple exclusion process (ASEP). Specifically, we con- 
sider the case p = 0.5 and L = 64. \ an d G are cal- 
culated theoretically as x{L) — 0.25 and <r(L) — 0.043. 
The diffusion constant D(L) is measured numerically as 
D{L) = 0.68 (see. Fig. 6 is also measured as 

®{L) = 3.9 from (Q2). Thus the relation D(L) X {L) = 
a(L)Q(L) seems valid. Study on the large scale behavior 
of the system might give new insights on ASEP. 

Third, we remark a connection with fluctuating hy- 
drodynamics. Note that the system in our numerical 
experiments may be too small for the dynamical behav- 
ior to be described within the framework of fluctuating 
hydrodynamics, while large enough so that we can check 
universal relations around equilibrium states. The fact 
that D(L) oc \f~L (See FigEJ m such small systems seems 
to imply that bare parameters appearing in fluctuating 
hydrodynamics depend on a length scale I at which hy- 
drodynamic description is defined. (Note that the scaling 
relation D(L) cx \f~L can be derived within the framework 
of fluctuating hydrodynamics, irrespective of values of 
the bare parameters |5j.) Thus, it is plausible to conjec- 
ture that we need to specify a scale I explicitly so as to 
have a hydrodynamic description. We do not understand 
how the formal derivation in Ref. 01 corresponds to this 
observation. 

Forth, we give a comment on diffusion of a marked 
particle. The diffusion is characterized by a tracer dif- 
fusion constant, which is not equal to the bulk diffusion 
constant, D, defined by (|10fl . One may has a question 
whether the tracer diffusion constant can be involved in 
a universal relation for NESS. We do not find such a re- 
lation yet, but it may be an interesting future problem to 
seek it. We note that the tracer diffusion constant does 
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not have the anomalous L dependence in contrast to the 
case of D. 

The final remark we wish to present regards the possi- 
bility of an experiment that can check the extended Ein- 
stein relations we have proposed. We conjecture that it 
may be best to study a system in which colloid particles, 
confined to one-dimension, are driven through a periodic 



array of optical traps [18| . By defining the effective tem- 
perature in the same way as did in (|12[) . the validity of 
the extended Einstein relations will be checked. 
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